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We present a theory for the estimation of a scalar or a vector magnetic field by its influence on an 
ensemble of trapped spin polarized atoms. The atoms interact off-resonantly with a continuous laser 
field, and the measurement of the polarization rotation of the probe light, induced by the dispersive 
atom-light coupling, leads to spin-squeezing of the atomic sample which enables an estimate of the 
magnetic field which is more precise than that expected from standard counting statistics. For 
polarized light and polarized atoms, a description of the non-classical components of the collective 
spin angular momentum for the atoms and the collective Stokes vectors of the light-field in terms 
of effective gaussian position and momentum variables is practically exact. The gaussian formalism 
describes the dynamics of the system very effectively and accounts explicitly for the back-action on 
the atoms due to measurement and for the estimate of the magnetic field. Multi-component magnetic 
fields are estimated by the measurement of suitably chosen atomic observables and precision and 
efficiency is gained by dividing the atomic gas in two or more samples which are entangled by the 
dispersive atom-light interaction. 



I. INTRODUCTION 

Precision atomic magnetometry relies on the measure- 
ment of the Larmor precession of a spin-polarized atomic 
sample in a magnetic field 0, 0, From standard 
counting statistics arguments, one might expect the un- 
certainty in such measurements to decrease with the in- 
teraction time t and with the number of atoms N at as 
l/y/Ne&t. If, on the other hand, the monitoring of the 
atomic sample, necessary for the read-out of the estimate 
of the magnetic field, squeezes the atomic spin, the above 
limit may be surpassed. In a recent theoretical analysis 
it was considered to estimate a scalar B field by a po- 
larization rotation measurement of an off-resonant light 
beam passing through a trapped cloud of spin- 1/2 atoms. 
This interaction squeezes the spin of the atomic sample, 
and by quantum trajectory theory combined with the 
classical theory of Kalman filters 0, la] , the uncertainty in 
the field strength was found to decrease as l/(N at t 3 / 2 ) Q. 
Very recently this proposal was implemented experimen- 
tally, and indeed sub-shotnoise sensitivity was found @. 

In a recent analysis of the experiment, we advocated 
treating all variables, including the magnetic field, as 
quantum variables @|. Secondly, we introduced a gaus- 
sian approximation at an early stage in the formula- 
tion of the theory. We further motivated and devel- 
oped this point of view in a detailed discussion of the 
spin-squeezing process @. As mentioned in these works, 
the gaussian approximation is essentially exact for the 
atomic and photonic degrees of freedom of the system un- 
der concern, and the advantages obtained by introducing 
this description from the outset of the theoretical treat- 
ment are at least four-fold: (i) the gaussian description 
explicitly accounts for the dynamics of the system and 
its behavior under measurements through update formu- 
lae for the expectation values and the covariance matrix 
which together fully characterize the gaussian state, (ii) 
the numerical treatment of the update formulae involves 



only the manipulation of low-dimensional matrices, (iii) 
in the limit of small time-steps the update formula for 
the covariance matrix translates into a matrix Ricatti 
differential equation which often lends itself to analyti- 
cal solution, and (iv) effects of noise introduced by, e.g., 
photon absorption and atomic decay are readily included. 

Here, we extend our previous analysis 0] to explore the 
possibilities for estimating B fields with not only one, but 
also two or three spatial components. In cases with more 
than one component, it is advantageous to use two or 
more polarized atomic samples. With such setups, we 
may identify sets of commuting observables which allow 
a simultaneous estimate of the B field components. We 
also discuss how to gain precision and efficiency by en- 
tangling the atomic gasses. 

The paper is organized as follows. In Sec. |H] we in- 
troduce the atom-photon system used for the estimation 
of the magnetic field and describe the atomic and pho- 
tonic gaussian variables. In Sec. IIIII we investigate the 
estimation of a single B field component, we describe 
our theoretical method in some detail and we derive an 
analytical solution for the decrease in variance of the B 
field as a function of time. In Sec. IIVI we present our 
results for the estimation of two or three spatial B field 
components. In Sec. we quantify the entanglement 
between the samples used in our optimal protocol for the 
estimation of several B field components. In Sec. lVH we 
explain how to include noise in the description and we 
study the effects of noise on the precision of measure- 
ments. In Sec. lVIIl we conclude and present an outlook. 



II. ATOM-LIGHT SYSTEM: COLLECTIVE 
VARIABLES 

To estimate the strength of a B field, we let it interact 
with an atomic spin-system which is continuo usly probed 
by a light beam along the lines of Refs. HH El El El 
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ILH, H3 |. In short, we imagine to have a gas of trapped 
spin- 1/2 atoms which are described by a collective spin 
operator J = § J2i 17 i wr th cr, the Pauli spin matrices. 
The atoms are initially pumped such that they are po- 
larized along the x axis and J x can be treated as a classi- 
cal variable (J x ) — -fy*- with N at the number of atoms. 
The other two projections of the spin, J y and J z , obey 
the commutation relation [J y ,J z ] = ihJ x which may be 
rewritten as [x a t,p a t] = * f° r the effective position and 



momentum variables x a 



= iPat 



= . The 



uncertainty is easily shown to be minimal in the initial 
state and, hence, the state pertaining to x a t and p a t is 
gaussian. 

The light beam propagates along the y axis and is 
linearly polarized along x such that its Stokes operator 
(S x ) = rt ^ ph is classical with N p h the number of photons. 
The two remaining components fulfill a commutator re- 
lation similar to the atomic spin case. Accordingly, for 
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the effective variables x D h = ; v , p P h = , , we 

have [x p h,Pph] = i and the initial coherent state of the 
field is a minimum uncertainty gaussian state. 

As shown in Ref. [3] and references therein, the light 
and atomic variables evolve as 



X a t l— * x at + K rPphi 
Xph > K rPat + X p h, 



Pat 
Pph > 



"Pat 
Pph 



(1) 

(2) 



when a segment of probe light of duration r and flux 
<&, and w ith a characteristic atom-light coupling k t oc 
a/ (J x )&t, is transmitted through the gas. The present 
dynamics in combination with a detection of the x p h com- 
ponent of the photon field leads to squeezing of the p at 
component of the atomic spin. 



III. ESTIMATING ONE SPATIAL 
COMPONENT OF A MAGNETIC FIELD 

The problem of estimating a single component of a B 
field was treated in Ref. p| and only a brief discussion 
is included here for completeness. Figure shows the 
setup. The B field directed along the y direction, causes 
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FIG. 1: Setup for measuring the B field component along the 
y axis. This is done using one atomic gas polarized along the 
x axis and one photon probe beam propagating along the y 
axis with a classical Sx- 

a Larmor rotation of the atomic spin towards the z axis, 
i.e., in time r, p at evolves as p at i— > p at — [i r B where \x T is 

given by the magnetic moment /?, via /i r = ^jt*-T @- 



Hence, Eqs. (!}-©) generalize to 
y S T y 

with y = {B y ,x a ,t,P & t,Xp] 1 ,Pp\ 1 ) T and 
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It is the coupling of the B field to the spin-squeezed vari- 
able Pat that makes an improved precision measurement 
of the magnetic field possible Q . 

Gaussian variables have been studied widely in relation 
to entanglement The gaussian description can be 

used as long as the interactions in the system are at most 
second order polynomials in the position and momentum 
operators and only homodyne measurements are carried 
out on the observables In particular, the operations 
of relevance for this work preserve the character of a gaus- 
sian state. We treat the continuous probe beam as a suc- 
cession of small beam segments of duration r, and let the 
state develop during time r between successive measure- 
ments. We recall that a gaussian state is fully character- 
ized by its mean value vector m = (y) and its covariance 
matrix 7 where jij = 2Re((y; - (yi))(yj - (yj)))- Accord- 
ingly, we only need update formulae for m and 7. The 
initial covariance matrix is 70 = diag(2 Var(_Bo), 1, 1, 1, 1) 
with Var(i?o) the initial variance of the B field. Under 
the linear transformation QJ, m and 7 transform as 



m(i + t) = S r m(<) 
7(£ + T) = S r7 (i)S^. 



(5) 
(6) 



The photon field is monitored continuously by detec- 
tion of x at . A major advantage of the gaussian descrip- 
tion is that the back-action on the residual system due to 
measurement is explicitly given. We write the covariance 
matrix as H El 111 



A 7 C 7 



(7) 



with A 7 the covariance matrix for the B field and atoms, 
yi = (By, Xnt,Pa,t) T , B 7 the covariance matrix for the 
photons, y2 = (x p h,Pph) T , and C 7 the correlation matrix 
for yi and yj . The measurement of x p h then transforms 
these matrices according to [TtI IT3 | 



A 7 — C 7 (7rB 7 7r) C 



B, 



*-7 

12x2 



0. 



(8b) 
(8c) 



where it — diag(l, 0), and ()~ denotes the Moore-Penrose 
pseudo-inverse. Equations l|8b|l and l|8cji follow from the 
fact that a new light segment is used in every measure- 
ment: The initial covariance matrix for every new seg- 
ment of a coherent photon beam is the identity matrix, 
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the resulting solution for the variance of the B field reads: 



FIG. 2: Uncertainty of one B field component as a function of 
time. The value at t = 5 ms is AB y = 5.814 x 10~ 5 pT. We 
have chosen a segment duration r = 10~ 8 s and corresponding 
field parameters = 0.0183 and fi T = 8.8 x 10~ 4 . 



and immediately after measurement there are no corre- 
lations between yi and y2 so C 7 is substituted with the 
zero matrix. 

The time evolution of m depends on the actual mea- 
surements in the optical detection which is a random pro- 
cess. Therefore the time evolution of the mean value vec- 
tor is a stochastic process, and it transforms as [^ ITfjlTij l 



mi i-> rri! + C 7 (7rB 7 7r) (x, ■ f 



(9) 



where \ is the difference between the measurement out- 
come and the expectation value of x p h, i.e., a gaussian 
random variable with mean value zero and variance V 2 - 
Since (7rB 7 7r)~ = diag(i? 11 1 , 0) with B\\ twice the vari- 
ance of a; p h, the second entrance in the vector (x, • ) need 
not be specified. In closing this section, we note that it is 
possible to understand the transformations in Eqs. jSajl 
and © by the corresponding transformation of classical 
gaussian probability distributions [9, 19]. 
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where we have introduced the photon flux $. 

FigureElshows the decrease in the uncertainty of the B 
field with time in a calculation with physically realizable 
parameters. From standard counting statistics one might 
have expected the variance to decrease as 1 /N at . Due to 
the measurement-induced squeezing of the atomic spin, 
we do, however, obtain the faster 1 /n^ ± decrease. As we 
shall discuss in Sec. IVII the inclusion of noise due to 
decoherence of atomic spins will alter this dependence 
on time and on the number of atoms. 



IV. ESTIMATING TWO OR THREE SPATIAL 
COMPONENTS OF A MAGNETIC FIELD 

In this section, we describe how to estimate two or 
three spatial components of a B field. In the setup in 
Fig. H we obtained an estimate of one component by 
using one atomic gas and one probe beam. To estimate 
more components we shall need more probe beams and 
more atomic gasses. 

In order to make a fair comparison of different schemes, 
we shall assume that all measurements are carried out 
in a time interval of the same duration, e.g., 5 ms as in 
Fig. 13 and that the total photon number used and the 
total number of atoms are kept constant. 



Two components: two probe beams and one or 
two separate atomic gasses 



A. Analytical solution 

From the update formulae, we may, in the limit of small 
time increments, derive a differential equation for A 7 . As 
By only causes rotation perpendicular to its direction, 
the variable x a t oc J y does not couple to (B y ,pat) and, 
hence, we only need to consider a 2 x 2 system with y = 
(By : p at ) T . The pertaining differential equation is on the 
matrix Ricatti form [2f| 

A 7 (t) = C- DA 7 (t) - A 7 (i)E - A 7 (i)BA 7 (t) , (10) 

with C = 0, D = (00), E = D r , and B = (° k ° 2 ) 
where n 2 = k^/t and /i = /i T /r. As may be checked 
by insertion, the solution to Eq. 11 Oil is A 7 = WU _1 , 
where W = -DW + CU and U = BW + EU. The 
linear differential equations for W, U can be solved, and 



In order to estimate B z in addition to B y , we observe 
that B z causes a rotation of x&t oc J y so if we add a second 
optical probe beam propagating along z in Fig. ^ an d 
with S x classical then x a t will cause a rotation of the field 
variable p p h oc S z on the second beam which we can then 
measure. The setup is symmetric with respect to y and 
z, and we obtain equal uncertainties on B y and B z . The 
problem with this approach is that unlike the unknown 
classical quantities B y and B z , the atomic observables x at 
and pat do not commute. While the first beam squeezes 
x a t, it anti-squeezes p a t, and the other beam does the 
opposite. Thus effectively we have no squeezing of the 
atoms leaving us with the 1 /t decrease in the variance of 
the B components as shown in Fig. 

To estimate simultaneously and precisely two B field 
components we need to measure two commuting atomic 
variables. Such a measurement is possible by using two 
separate gasses and by estimating one B field component 
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FIG. 3: Uncertainty of two B field components as a function of 
time. We use a single atomic gas and two probe beams which 
are turned on simultaneously. We have chosen a segment 
duration r = 1CP 8 s and corresponding field parameters — 
0.0183 and fJ, T = 8.8 x 10~ 4 . The joint uncertainty of the two 
B fields at t = 5 ms is AB y = AB Z = 0.1083 pT. 



on each system with individual probe beams. Both sys- 
tems are then equivalent to the setup in Fig. used to 
measure B y , but in one system (not shown), the probe 
beam propagates along the z axis such that we estimate 
B z by measuring p p h- The result is similar to the result 
in Fig.|2except that the number of atoms and the photon 
flux are both divided by two as they are shared between 
the two systems. The uncertainty of the B fields for large 
t is proportional to W^ft** 3 so the uncertainty for esti- 
mating two B field components is 2\/2 times larger than 
if we had measured only one with the same field and 
atomic resources. 

The achievements of a sequential measurement of first 
By and then B z are shown in Fig. 0J In the first half of 
the time the full photon flux is spent to measure B y , leav- 
ing the uncertainty about B z unchanged. Although B z is 
subsequently coupled to an anti-squeezed atomic compo- 
nent, it is quickly squeezed, and the estimate of B z shows 
a W^at** 3 dependence as in Eq. ifTTTl . Now the other 
atomic observable is anti-squeezed, but this will surely 
not degrade our information already obtained about the 
classical B y component. The values of Vax(B y (t)) and 
Va,r(B z (t)) are 2\/2 times larger than the uncertainty re- 
ported in Fig. [21 because only half the time is spent on 
the measurement of each component. 



B. Two components: two entangled gasses and two 
probe beams 

If we split the atomic sample into two and polarize one 
gas along x and the other along —x such that (J Xl ) — 
— (Jx 2 )i then the two observables (J Vl — J y2 ) and (J Zl — 
J Z2 ), and equivalently x atl -x a t 2 andp atl -p a t 2 commute. 
We couple these observables to the B fields and probe 
beams and use the setup shown in Fig. (see also the 
entanglement experiment of Ref. [T^l. 
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FIG. 4: Uncertainty of two B field components as a function 
of time using one atomic gas. B y is estimated in the first half 
of the time and B z in the second half. The dashed line is 
for B y and the value at t = 5 ms is AB y = 1.647 x 10~ 4 pT. 
The full line is for B z and the value at t = 5 ms is AB Z = 
1.645 x 10 _4 pT. We have chosen a segment duration r = 
10 -8 s and corresponding field parameters = 0.0183 and 
fj, T = 8.8 x 10 -4 . 



Both optical probe beams have S x classical, and one 
beam propagates along y, the other along z. The beams 




FIG. 5: Setup for measuring two components of a B field 
using two entangled gasses and two probe beams. 

pass through both gasses and, oppositely to the protocol 
in the previous section, we use all atoms to estimate the 
B field components. The effective Hamiltonian for the 
setup is 



HmtT = H T B y (x atl + X at2 ) + fl T B z (p atl + p a t 2 ) 

+ KriPa-U -PatJPphj + K T (x atl - Xat 2 )£ph 2 , 

(12) 



where the two minus signs can be implemented by chang- 
ing the sign on n T after the probe beams have passed 
through the first gas. In practice the change in sign 
can be effectuated by changing the sign of the detun- 
ing or by interchanging tx + and c~ polarizations with, 
e.g., a half- wave-plate [2l|. The gaussian state vector is 

y {-B z: B y , X a \ >1 , Pati i *£at 2 i J*at 2 ) '^ph 1 ! Pph 1 : *^ph 2 ' Pph 2 ) 

and from the Heisenberg equations of motion for the op- 



erators we get the following transformation matrix 
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The time evolution of the uncertainty of the B fields is 
shown in Fig. The final uncertainty of B y and B z is a 
factor \[2 higher than in Fig. |2 where we used the entire 
photon flux to probe only a single B field component, 
but a factor of two lower than in the setup with separate 
probing of non-entangled gasses. These factors were ex- 
pected because of the 1 /y/N% t <s>t 3 dependence of Eq. (11 II . 
By using entangled gasses we use all atoms, but only half 
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FIG. 7: Setup used to obtain estimates of all three compo- 
nents of a magnetic field. We use six atomic gasses and three 
probe beams. 
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FIG. 6: Uncertainty of two B field components as a func- 
tion of time using two atomic entangled gasses and two probe 
beams. We have chosen a segment duration r = 10 -8 s 
and corresponding field parameters k z t = 0.0183 and fj, T = 
8.8 x 10" 4 . The uncertainty at t = 5ms is AB y = AB Z = 
8.221 x 10~ 5 pT. 

of the photon flux to estimate each B field component. 



C. Three dimensional vector magnetometry 

For the estimation of three components of a magnetic 
field, the situation changes since a spin polarized sample 
is not a probe for the field component parallel with the 
spin. For the sequential probing one would thus measure 
By and B z as just described, but one would have to rotate 
the sample by 90° to determine the last component B x , 
and the errors of such a rotation will limit the precision. 

As an alternative, we can divide the gas into three 
samples which are spin polarized along different direc- 
tions. With such a setup, we can perform independent 
estimates of the three components. Since each compo- 
nent is determined by one third of the atoms and one 



third of the photons the scaling of Eq. Ill 111 predicts that 
the uncertainty is a factor of 3V3 larger than in Fig. [3 
where we only estimated one B field component. 

We wish to couple the B field components to three 
commuting atomic operators, involving as many atoms as 
possible, and a better, but also more complicated, setup 
is shown in Fig. Here we use six entangled gasses 
and we let three probe beams pass through four gasses 
each. The gasses are polarized with the following macro- 
scopic components: (J Xl ) = —(Jxi), (Jy 2 ) = ~{Jy$)> an d 
(Jzs) — ~(Jze)- The optical fields are linearly polar- 
ized with macroscopic Stokes parameters (S Z1 ), (S V2 ), 
and (S X3 ). The Hamiltonian for the system shown in 
Fig. is given by 



HmtTy/h(J Xl ) = {J, T (Jx 2 + Jx 3 + Jxs + Jx 6 )B x 

+ ^ri'Jyi + Jy 3 + Jy 4 + Jys)B y 
+ Ht(Jzi + Jz 2 + J Zi + Jz a )B z 

+ K T (J Z2 — Jy 3 — J Z& + Jy 6 )S^ 

+ K T (J Z1 ~ J X3 — J Z4 + J X(i )S2 
"b ^r(Jyi J x 2 <A/4 ~t~ Jx^^l: 

Terms like [i T J xl B x , which couple the classical com- 
ponents of the atomic spins to the B fields, are omit- 
ted from the interaction Hamiltonian as they do not 
contribute to the interactions to the same order, e.g., 
[J yi , H T Jx 1 B x ] — —i[i T hJ Zl B x , the product of two small 
quantities. By using the Hamiltonian of Eq. 11141 we mea- 
sure three commuting observables J Z2 — J V3 — J Z5 + J ye , 
J Zl - Jx 3 - J Zi + Jx 6 , and J yi ~ J X2 - J yi + J X5 and from 
their commutators with the Larmor term in Eq. 11141 , we 
see that they evolve in direct proportion with the three 
B field components. 

The uncertainty of the B field components is shown as 
a function of time in Fig. |S| As we use 4 /6 of the atoms 
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t = 5 ms are AB X = AB y = AB Z = 1.510 x 10" 4 pT. 



and !/3 of the photon flux to estimate each B field com- 
ponent, the value of the uncertainty at t = 5 ms is 3 ^/2 
times larger than if we use all atoms and all photons to 
estimate only one B field component and a factor of two 
smaller than if we use three separate systems to estimate 
the three B field components. 



QUANTIFYING ENTANGLEMENT 
BETWEEN ATOMIC SAMPLES 



To measure two or three B field components most effi- 
ciently, we showed that one should use entangled gasses. 
Here we quantify the degree of entanglement by calculat- 
ing the gaussian entanglement of formation (GEoF) j23| 
from the covariance matrix of the atoms 7 at0 ms- Every 
time we apply the update formula, we may extract 7 atoms 
for a pair of gasses from our numerical procedure and up 
to local unitary operations this matrix turns out to be 
on the form 



7a 



(n k x 

n -fcj 

k x n 

\0 -kp n 



(15) 



where n, k x , and k p (k x — k p , in our case) are the 
quantities of interest for the evaluation of the GEoF, 
E = c+(A) log[c+(A)] -c_(A) log[c_(A)], with c±(A) = 
±(A-V 2 ± A 1 / 2 ) 2 , and A = min (l, y/(n - k x )(n - fc p )) . 

The GEoF for the two gasses used to estimate two 
B field components is shown in Fig. El For three B 
field components, we used six atomic gasses, and we have 
calculated the GEoF between different pairs of the gasses. 
Figure El shows the GEoF between two gasses polarized 
in opposite directions, e.g., gas number 1 and 4 in Fig. [HI 
The GEoF between pairs like 1 and 2 is zero. 

The setup with two gasses is quite equivalent to the one 
implemented in recent entanglement experiments |l2l | ex- 
cept that the atomic systems are under the additional 




FIG. 10: GEoF for two gasses polarized in opposite directions 
as considered in Fig. El 



influence of an initially unknown B field. This slows 
down the initial rate of generation of entanglement, but 
as Var(S(t)) approaches zero, the entanglement grows 
without limits as long as absorption and atomic decay 
can be neglected [131 ■ In the case of six gasses which 
are probed in a non-symmetric way, some pairs show en- 
tanglement and some do not. This can be understood by 
identification of operators that do not couple to the probe 
fields. The convergence of the entanglement between op- 
positely polarized gasses towards a constant value is also 
observed without coupling to a B field, and it is due 
to the incompleteness of the measurements on the pair 
by fields that also couple to other pairs of gasses. In 
symmetric setups with multiple gaussian variables, the 
theoretical maximum of pairwise entanglement between 
systems also have upper limits reflecting the impossibil- 
ity for a quantum system to be maximally entangled with 
several other quantum systems at the same time (24ll25| . 
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VI. EFFECTS OF NOISE 

Effects of noise were recently discussed in Ref. Q. 
Here we include noise in our gaussian description. As 
the photonic probe beams pass through the atomic sam- 
ple, there is a probability for photon absorption j2^| 

and a related probability of atomic 



N, 



7« 



74 



7«+a ; 



■, where T is the atomic 



at a r a /4+A 2 

decay |Rl2fil| 77- = 

decay rate, a — A 7( 27r ) lS the resonant photon absorp- 
tion cross-section, A is the beam cross-section, and A 
the detuning. These processes lead to a reduction in the 
polarization of the Stokes vector and the atomic spin and 
to an incoherent noise contribution. These features were 
discussed in Refs. 0,113] and at length in Ref. @ so here 
it is sufficient to recall the generalizations of the update 
formulae in Eqs. 



va{t + r) = L r S T m(i) 
7 (t + r) - L r S T7 (i)S^L r 



M, 



hN ph 



2(5 x (t)) 



(16) 
N, 



(17) 



where the diagonal matrices L r = 
diag( 1 , V 1 - Vt, V 1 - Vt, VI - e, VI - e) , M T = 
diag(0, rj T ,r] T , 0, 0) and N = diag(0, 0, 0, e, e) describe the 
loss of polarization, and the noise introduced by stimu- 
lated emission and photon absorption, respectively. In 
each time step r the atomic polarization (J x ) is reduced 
by the factor (1 — r/ T ) and k t is reduced by \J\ — rj T . 
The presence of noise leaves the update formula for mea- 
surements in Eq. © unchanged. We may analyze the 
pertaining Ricatti equation, which with y = (£?,p a t) T 
can be written on the same form as in Eq. llOli but 

now with C = (o_j4_,),D= (°5j),E = D T , and 

B = — f g) V? \ , where n = t) t /t. One can 

show that the noise terms in D and E will have vanishing 
effect, and if we restrict ourselves to times corresponding 
to r\t <C 1, we may neglect the time dependence of (J x ), 
(S x ), and k such that N °*/(J X ) = 2, Ar P h/<s x ) = 2, and 
K(t) = k(0). With these approximations, the resulting 
linear equations for the matrices U and W can be solved 
analytically, leading to lengthy expressions with sums of 
products of exponential functions, constant terms, and 
terms linear in time t. In the limit yj rjK 2 t ^ 1, it is 
an accurate approximation to maintain only the leading 
exponential and we find 



Var(B(t)) - -L 



ZL 



(18) 



Compared with the result in Eq. l(TT|l . we note that in 
the long time limit the uncertainty decreases as 1 /N at t, 
and not as Vivf t t 3 as in the noise-less case. 

Figure El shows (solid line) the time evolution of the 
uncertainty of the B field. For short times the analytical 
expression 11 111 disregarding the noise (dotted line) in 



the figure, is a good approximation, and for long times 
the numerical result follows the expression l|T8jl . shown 
as a dashed curve on the interval between 0.5 ms and 
5 ms. Similar results are obtained in the simulation of 
measurements of two and three B field components. 



% °- 01 - 

< 0.001- 



10" 




t [ms] 



FIG. 11: Uncertainty of one B field component as a function 
of time. The full line is the result of a full numerical calcu- 
lation and the value at t = 5 ms is AB y = 2.333 x 10~ 4 pT. 
The dotted line is without inclusion of noise and the value at 
t — 5 ms is AB y = 5.814 x 10 -5 pT. We have chosen a seg- 
ment duration r = 10 -8 s and corresponding field and noise 
parameters «? = 0.0183, fjb T = 8.8 x 10 -4 , r) T = 1.76 x 10" 8 
and e = 0.0281. The dashed curve shows the result of Eq. USt 
valid for t <gC 1 s, and sJr\K?t S> 1. 



VII. CONCLUSION AND OUTLOOK 

We have considered how to estimate a vector magnetic 
field using a gaussian description of the variables describ- 
ing the system. To estimate more than one field com- 
ponent it is fruitful to use pairwise entangled separate 
gasses. In other applications, e.g. teleportation, shared 
entanglement over some distance is a useful resource. We 
showed that entanglement can be a useful local resource 
to improve the accuracy in measurements and parameter 
estimation. We note that our protocol for the estimation 
of two B field components using two optical probe beams 
and two entangled gasses is experimentally very feasible. 
Essentially, it would require a combination of the mag- 
netometry setup of Ref. Q and the entanglement setup 
of Ref. E3- 

In a broader perspective, the present work brings out 
the virtues of the gaussian state formalism when it comes 
not only to the detailed characterization of entanglement 
(see, e.g., Refs. 0, 0, 0, 0] and references therein), 
but also to the practical description of quantum sys- 
tems 0, El, E3| ■ The fact that the gaussian state is fully 
characterized in terms of its expectation value vector and 
covariance matrix means that the theory is easy to for- 
mulate and to evaluate. It is an outstanding advantage 
of the gaussian state description that explicit update for- 
mulae exist not only for the interaction dynamics but also 
for the back-action due to measurement. The gaussian 
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state description is a very versatile tool and we foresee 
this approach to be used for the description of a variety 
of different systems in the near future. 



Acknowledgments 



L.B.M. is supported by the Danish Natural Science 
Research Council (Grant No. 21-03-0163). 



[1] D. Budker, W. Gawlik, D. F. Kimball, S. M. Rochester, 
V. V. Yashchuk, and A. Weis, Rev. Mod. Phys. 74, 1153 
(2002). 

[2] I. K. Kominis, T. W. Kornack, J. C. Alfred, and M. V. 
Romalis, Nature (London) 422, 596 (2003). 

[3] M. Auzinsh, D. Budker, D. F. Kimball, S. M. Rochester, 
J. E. Stalnaker, A. O. Sushkov, and V. V. Yashchuk, 
Can a quantum nondemolition measurement improve the 
sensitivity of an atomic magnetometer?, physics/0403097 
(2004). 

[4] H. Carmichael, An Open Systems Approach to Quantum 

Optics (Springer- Verlag, Berlin Heidelberg, 1993). 
[5] J. M. Geremia, J. K. Stockton, A. C. Doherty, and 

H. Mabuchi, Phys. Rev. Lett. 91, 250801 (2003). 
[6] J. K. Stockton, J. M. Geremia, A. C. Doherty, and 

H. Mabuchi, Phys. Rev. A 69, 032109 (2004). 
[7] J. M. Geremia, J. K. Stockton, and H. Mabuchi, 

Sub-shotnoise atomic magnetometry, quant-ph/0401107 

(2004). 

[8] K. M0lmer and L. B. Madsen, Estimation of a classi- 
cal parameter with gaussian probes: Magnetometry with 
collective atomic spins, quant-ph/0402158 (2004), Phys. 
Rev. A, at print. 
[9] L. B. Madsen and K. M0lmer, Spin squeezing and pre- 
cision probing with light and samples of atoms in the 
gaussian approximation, quant-ph/0406146 (2004), Phys. 
Rev. A, at print. 

[10] A. Kuzmich, N. P. Bigelow, and L. Mandel, Europhys. 
Lett. 42, 481 (1998). 

[11] L.-M. Duan, J. I. Cirac, P. Zoller, and E. S. Polzik, Phys. 
Rev. Lett. 85, 5643 (2000). 

[12] B. Julsgaard, A. Kozhekin, and E. S. Polzik, Nature 
(London) 413, 400 (2001). 

[13] Y. Takahashi, K. Honda, N. Tanaka, K. Toyoda, 



K. Ishikawa, and T. Yabuzaki, Phys. Rev. A 60, 4974 
(1999). 

[14] A. Kuzmich, L. Mandel, J. Janis, Y. E. Young, R. Ejnis- 
man, and N. P. Bigelow, Phys. Rev. A 60, 2346 (1999). 

[15] G. Giedke, J. Eisert, J. I. Cirac, and M. B. Plenio, Entan- 
glement transformations of pure gaussian states, quant- 
ph/0301038 (2003). 

[16] G. Giedke and J. I. Cirac, Phys. Rev. A 66, 032316 

(2002) . 

[17] J. Fiurasek, Phys. Rev. Lett. 89, 137904 (2002). 

[18] J. Eisert and M. B. Plenio, Introduction to the basics 
of entanglement theory in continuous-variable systems, 
quant-ph/0312071 (2003). 

[19] P. S. Maybeck, Stochastic Models, Estimation and Con- 
trol, vol. 1 of Mathematics in Science and Engineering, 
Volume 141 (Academic Press, New York, 1979). 

[20] J. K. Stockton, J. M. Geremia, A. C. Doherty, and 
H. Mabuchi, Robust quantum parameter estimation: Co- 
herent magnetometry with feedback, quant-ph/0309101 

(2003) . 

[21] W. A. Shurcliff, Polarized Light (Harvard University 
Press, Cambridge, 1962). 

[22] G. Giedke, M. M. Wolf, O. Kriiger, R. F. Werner, and 
J. I. Cirac, Phys. Rev. Lett. 91, 107901 (2003). 

[23] J. Sherson and K. M0lmer, In preparation (2004). 

[24] M. M. Wolf, F. Verstraete, and J. I. Cirac, Phys. Rev. 
Lett. 92, 087903 (2004). 

[25] M. B. Plenio, J. Eisert, J. Dreifiig, and M. Cramer, En- 
tropy, entanglement, and area: analytical results for har- 
monic lattice systems, quant-ph/0405142 (2004). 

[26] K. Hammerer, K. M0lmer, E. S. Polzik, and 
J. I. Cirac, Light-matter quantum interface, quant- 
ph/0312156 (2003), Phys. Rev. A, at print. 



